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INTRODUCTION 


This IS the final report covering the technical work per- 
formed by the Department of Aerospace Engineering and Applied 
Mechanics under Contract No. NAS2-7850 (U-S. Air Mobility 
Research and Development Laboratory - Moffett Field, California) . 
The work was supervised by Dr. Widen Tabakoff, Professor in 
the Department of Aerospace Engineering and Applied Mechanics. 

The report covers work conducted from October 1, 1973 to 
October 30, 1977. 

The contents of the technical reports written during the 
contract period have been incorporated in this report through 
the provision of abstracts. In addition, the technical 
publications in Scientific Journals and Papers are given with 
the first page of the corresponding publication. During the 
last year of this contract we started to develop a method 
of analyzing turbulent nonadiabatic viscous flov; through radial 
turbomachine rorors. This work is now only in the progress 
stage . 
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Abstract 


A numerxcai study to determxne the 
temperature distribution in the rotor of a 
radial inflow turbine is presented- The 
study IS based on the use of the finite 
element metnod in the three dimensional 
heat conduction proolem. Different cool- 
ing techniques witn various coolant: to 
primary mass flow ratios are investigated. 
The resulting temperature distribution in 
the rotor are presented for comparison. 
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University of CincmnatL 
C^ncmnatL Ohio 


losses in Radial InHow Mines 

A dudy wa$ conducUd to determine erpenTneTiially and theoretically the losses 
radial inflow turbine 71022 ^ 5 . Ezienswe experimental dcUa ivere obtained to investigate 
the flow behavior in a fulUscale radial turbine stator annidus, *4 theoretical model to 
predict the losses in both the vaned and vaneless regions of the nozzle was developed. In 
this analysis, the interaction effects between the stator and the rotor are not considered. 
It was found that the losses incurred due to the end wall boundary layers can he cignifi^ 
cant, especitdly if they are characterized by a strong crossflow. The losses estimated 
vsnng the analytical study are compared with the experimentally determined values. 


Introduction 

Due to the complesd-cy of the Slow patterns m radial inflow 
turbines, previous analytical investigations dealt with ideal flow, 
with relatively simple loss models. The viscous losses were as- 
sumed to be proportional to the average kmetic energy m refer- 
ence [1],^ while approxmiatioas based on flat plate and pipe 
flow analogies were used m reference [2] 

In the present study, experimental measurements were ob- 
tamed for the flow at the inlet and the exit of a full-scale radial 
turbine stator annulus Detailed data, which includes flow direc- 
tions, total and static pressure surveys, have been used to gam 
an insight mto the flow processes by which the dostnbution of the 
losses is aflected. Results are then used to develop a better 
theoretical model for loss prediction m both the vaned and the 
vaneless neld constitutmg the stator annulus. It is hoped that 
this model will assist in developmg an optmuzatiou scheme which 
can be used for designing the stator annulus of a radial turbine 

Experimental Investigation 

A test facility was designed to mcorporate a complete radial 
inflow turbine for stator and rotor testing usmg cold air The 
present phase of expenmental mvestigation was concerned with 
evaiuatmg the overall performance of the stator annulus and to 
gam a better understanding of tbe flow behavior m the radial 
nozzle This was accomplished by testmg a fuU^cale stator as- 
sembly installed m the scroll housmg of the turbme. The rotor 
was not assembled during these tests. Instead, a plastic body of 
revolution was installed to provide a smooth, contmuous hub 
flow-path downstream the stator. The scroll which was located 
upstream the stator assembly was of the nonsymmetnc circular 
type. The nozzle stator annulus consisted of 13 untwisted con- 
stant section vanes of 0 04S2-m chord length The annulus had an 



^Numbers in brackets designate Hefereaces at end of paper 

Contributed by the Fluids Engineering Division of Ths Ambbic Society 
OF Msca.iMCAL Evgimiebs and presented &z the Gas Turbme and Fluids 
Engineenag Conrerence, New Orleans, La . Mar^b 21-25. 1&76 Macuscnpt 
received at A5ME Headquarters January 5. 1976 Paper No 76-FE-9 


outer diameter of 0 191 m and an inner diameter of 0 139 m 
The distance between the two side walls was constant and equal 
to 0.0147 m. The facility, instrumentation, testing techniques 
and the scator design characteristics are described in reference 
(3] Fig 1 shows a oross-sectional sketch of the stator assembly 
Surveys of total pressures, static pressures and the flow direc- 
tions were made both upstream and downstream of the stator 
annulus, using calibrated wedge probes Outside the boundary 
layer regions, a three-dimensional cobra-probe was designed bo 
measure the flow total pressure, static pressure, pitch and yaw 
angles The probe consisted of flve 0 037 cm outside diameter 
tube with their meaunng ends grounded to a total included 
angle of 70'* Inside the boundary layer regions a pitOD-static 
yawmeter was used to measure the air velocity and flow angle? 
The yawmeter tube tip was flattened to approximately 0 Glo 
cm inside height Measurements of vane surface pressure dt^- 
tnbutions were also obtained by means of ports of 0 045 cm 
diameter drilled m h 3 podernuc tubes embedded in the body of 
two consecutive vanes The upstream measuring station, 0, of 
Fig 1 was located 1/4 chord radially outward from the leading 
edge, while the downstream station 2 was at 0 1 chord radially 
inward from the trailing" edge (very close to station 1) Survey^ 
were obtained through the spanwise and picchwise movements 
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Heat Transfer in Cooled Guide Vane Blades 

W. Tabakoff, R. Kotwal and A. Hamed, 

University of Cincinnati, Cincinnati, Ohio 


Abstract 


A numerical study x:o determine the 
temperature distribution in the guide 
vanes of a radial inflow turbine is pre- 
sented. A computer program nas been de- 
veloped to calculate the temperature dis- 
tribution when the vanes are cooled inter- 
nally using a comoination of impingement 
and film cooling techniques. The study is 
based on tne use of the finite difference 
method in a two dimensional heat conduc- 
tion proolem. The results are tnen com- 
pared ho deternine the best cooling con- 
figuration for a certain coolant to pri- 
mary mass f lov/ ratio , 



Los Angeles, Calif /January 24-26, 1977 


For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics^ 
1290 Avenue of the Americas. New York, N Y 10019 



Reprinted from JOURNAL OF AIRCRAFT, VoL 14, No 2, February 1977, 173- L76 
Copyright, 1977, by the Amencan Institute of Aeronautics and Astronautics, and reprinted by permission of the copyright owner 


Temperature Distributioa Study in a Cooled 
Radial Inflow Turbine Rotor 






A. Hamed,* E. Baskharone.t and W. TabakoffJ 

University of Cincinnati, Cincinnati, Ohio 




& 


In this work, the three dlmenstonai temperature distnbution m the cooled rotor of a radial inflow turbine is 
determined numencally using the finite element method. Some of the results obiaiaed using different cooling 
techniques with vanous coolant to primary mass flow ratios are presented for companson. In this analysis, the 
complicated geometnes of the hot rotor and coolant passage surfaces are handled easily, and the temperatures 
are determined without loss of accuracy at these convective boundaries. The present work can easily be used tn 
combination with a finite element stress analysis, to investigate the thermal stresses corresponding to the dif- 
ferent cooling arrangements. 




C/ = 
D • = 
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Q 
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Re 

[SR] = 
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r 

V 
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V =S 

0 ) - 

Subscripts 
X,y,Z =: 


. ^ Nomenclature 

' skm fnctron coefficient 
internal cooling passage hydraulic diameter (ft) 
heat gener^bn rate per unit volume (Btu/ft^ hr) 
heat generation matrix, Eq. (9) 
convective heat transfer coefficient 
(Btu/ft hr*F) 
surface convection matrix 
coefficient of thermal conductivity (Btu/ft hr "iO 
thermal stiffness matrix 
Nusselt number 

outward normal unit vector from the rotor surface 
Prandtl number 

row vector representing the spatial distribution 
within the element, Eq. (5) 
heat flux density (Btu/ft^ hr) 
element nodal spatial matrix, Eq. (6) 

Reynolds number 

thermal boundary load matrix, Eq, (9) 

rotor surface exchanging heat with the hot gas or 

coolant flow 

rotor surfaces on which heat flux is specified, in- 
cluding adiabatic surfaces 
temperature 

the volume of the tetrahedral element m the finite 

element representation (ft^) 

column vector of coefficients, Eq, (5) 

kinematic viscosity (ft^/hr) 

rotor angular velocity (radiance/hr) 

refer to panial denvation m these directions, 
respectively 

flow conditions generally m hot gas or coolant 
flow 


Superscripts 
t 


= transpose of a tensor 


Introduction 

I n the last decade internal cooling has been used extensively 
resulting m the continuous increase of turbine inlet 
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temperatures Most internal cooling systems are designed 
using senu-empiricai methods to achieve the highest possible 
effectiveness. Reference 1 presents a review of the present 
state-of-the-art for the internal cooling of turbine nozzles m 
aircraft applications. 

Expenmentai studies of air cooled radial inflow turbine 
rotors are reported in Refs. 2 and 3 Eranger’ investigated the 
effeaiveness of veil cooling on the hub side of the rotor. It 
was found that the cooling effectiveness was larger at the 
rotor tip, and decreased as the cooling film is heated and 
mixed with the hot turbine flow Petrick and Smith ^ studied 
the effect of cooling the backside of the rotor on us tem- 
perature distnbution. While veil and rotor backside cooling 
techmques are effective for the rotor disk, they induce little 
vanation in the blade temperature. Internal cooling on the 
other hand can be used to produce the desired reduction in the 
blade temperature as reported m Ref. 4 

Considering the time and cost associated with expenmentai 
studies, computational methods can provide valuable design 
and development information, concermng the temperature 
and stress distnoutions associated with the different turbine 
cooling techniques. Several studies are available for the tem- 
perature distnbution in the cooled blades of the axial gas flow 
turbine blades. They are generally based on two dimensional 
temperature computations with a spanwise vanation m the 
blade surface convective heat transfer coefficient. In the 
radial inflow turbine rotor, a similar approacn is not possible 
and the real three dimensional temperature field must be in- 
vestigated. The finite element method is particularly suitable 
for handling such complicated geometries without loss of ac- 
curacy at the boundaries. 

In the present study a numerical method is developed to 
determine the temperature field in the rotor of a cooled radial 
inflow turbine The computations are based on the use of the 
finite element method, with a variational statement of the 
three dimensional conduction problem. The vanous boundary 
conditions, at the different external and internal cooling sur- 
faces of the rotor are discussed. Different cooling methods 
and coolant flow rates are investigated and the resulting 
cooling effectiveness and temperature distributions are re- 
ported, The data obtained from the present analysis were 
found to be m agreement with the available experimental 
.measurements. 


Governing Equations 

The field equation for the steady state three dimensional 
heat conduction problem in an isotropic medium can be ex- 
pressed as 
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With increased turbine inlet teippera- 
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stresses associated with different 
cooling patterns in a. radial inflow rur-- 
bine rotor* The finite element method 
IS used in the stress calculations taJcing 
into consideration centrifugal, thermal 
and aerodynamic loading. The effects 
of temperature distribution and the 
presence of internal cooling passages 
are discussed. 
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VISCOUS PLOW ANALYSIS IN TURBOi'4ACHINE ROTORS 


I.M. Khalil and W. Tabakoff 

A method of analyzing the turbulent nonadiabatic viscous flow 
through txirbomachine rouors is now in the progress stage. The 
field analysis is based upon the numerical integration of the full 
Navier-Stokes equations/ together with the energy equation, over 
the rotor b lade ~to-b lade stream surfaces. The numerical code 
used to solve the elliptic governing equations will employ a 
curvilinear body-oriented coordinate system that suits the most 
complicated blade geometries. Turbulence will be represented by 
a two-equation model, one expressing the development of the turbulent 
kinetic energy, and the other its dissipation rate. 

The flow analysis will be applied to various types of 
turbomachine rotors. The flow characteristics within the 
different rotor passages will be investigated over a wide range 
of operating conditions. The development of the viscous boundary 
layers over the blade surfaces will be also studied for different 
boundary conditions. Calculated results will be compared to 
experimental data. It is felt that the method of this analysis 
is quite general and can deal with wide range of applications. 

‘Preliminary computations, using the computer programs 
developed during this period of study, were performed for 
different flov/ configurations. For example, the viscous flows 
within the rotors of a radial bladed compressor and a mixed 
flow turbine are investigated under different operating conditions. 
Excellent results were obtained when compared with existing 
experimental data over a wide range of Reynolds numbers. 

Currently, active steps are being carried out to incorporate 
Thompson approach (Ref. 1) for transformation in our main code 
such that the complicated blade geometries can be dealt with. 

Up till now, we are faced with minor problems in the transforma- 
tion process and procedures are being devised to solve these 
problems. Specifically, the mesh resolution near the leading 
and trailing edge needs special attention. The Jacobian of 
transformation for these regions varies largely and leads to 
less accurate results if small step sizes are not considered. 



FL'OW MfMiYSIS 


The partial differential equations that govern the flow 
behavior within the turbomachine rotor passages are presented. 

A transformation of the three dimensional flow equations to 
several particular two dimensional forms, on predetermined 
stream surfaces, is then outlined. The details of the mathematical 
formulation used to perform the transformation is not reported 
here, but the final results are summarized. The resulting flow 
equations are further expressed in the form of conservation 
law by using the vorticity-stream function formulation. Finally, 
the boundary conditions necessary to obtain a unique solution 
ro the problem are discussed. 

Fundamental Aero thermodynamic Relations 

The three dimensional flow of viscous, compressible fluid 
through turbomachinery is governed by the following set of laws. 

Conservation of Mass ; 

|| + V- (pW) = 0 (1) 

Conservation of Momentum ; 

For a blade rov7 rotating at a constant angular velocity 
about its axis, Newton's second lav/ of motion when combined with 
Stokes hypothesis can be written as 

P (It + W-VW + 20xW - = - Vp - Vx[y (V X W) ] 

+ J V(p7*W) (2) 

Conservation of Energy ; 

In the absence of heat sources or sinks, the first law of 
thermodynamics for a heat conducting fluid with a thermal 
conductivity (K) can be v/ritten as : 

p (|~ + W- Vh) = ||- + (W- V) p + J) + V • (KVT) {3) 
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where 


p = static pressure. 


p = static density, 

T = statxc temperature, 
h = static enthalpy 
W - relative velocity vector, 

R = position vector 

V = kinematic viscosity coefficient, 
jp = viscous dissipation function. 


The relations among the flow state variables are those of the 
ideal gas equations ; 


p = pRT 

and dh = C dT 

P 


OF POOR 


( 4 ) 

( 5 ) 


In order that the above equations represent a turbulent flow 
situation, the mass-averaged variables principle is utilized. 
According to this principle, the components of the velocity 
vector W, and the temperature T, are expressed on a mass-averaged 
basis, while the pressure and the density p are expressed 
by their mean values only. For any variable QCx^,X 2/X2, tl , 
its mass-average is given as follows: 


q(x^,X2,x^) = Q(x^,X2,X2,-c) P (x^ , x^ , x^ , t) /p 

where P (x^,X2,X2,t) is the instantaneous value of the density at 

any point whose location is defined by the coordinate (x^,X2,X2). 

The over bar in the above equation represents the conventional 

time average, thus P = p. An effective viscosity u is also used 

0 

instead of the kinematic viscosity y . This effective viscosity 
IS assumed to describe the effecrs of Reynolds stress. In the 
present work, a two equation turbulence model will be employed 
to calculate 
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Preliminaries 


An approach is taken to reduce the spatial dimensions of the 
problem from the general three dimensional form to several 
particular two dimensional forms. The object is to obtain a 
solution to the flow governed by the equations (1) to (5) 
through an appropraate combination of two dimensional flow 

solutions. The reduction in spatial dimension is achieved 

\ 

through a consideration of the blade-to-blade stream surface 
concept. The blade-to-blade stream surface/ S^/ may be 
described as the surface that would extend from the pressure 
surface of a blade to the suction surface of the next blade 
as shown in Figure la. If desired, this surface may be considered 
to be a surface of revolution generated by a meridional stream- 
line obtained from a meridional flov/ analysis. On this srream 
surface, the flov; equations Cli , (2) and (3) are transformed to 

several two dimensional mathematical expressions. This results 
from the fact that the surface provides relations amongst 
the coordinare variables, such that the variation of flow 
properties over the surface may be described in terms of tv/o 
space variables only. If a solution to the resulting two 
dimensional equations is obtained, then v/ith several such blade- 
to-blade solutions, the flov/ properties throughout the three 
dimensional field may be evaluated. 

Stream Surfaces Equations ; 

Consider a flow annulus as shown in Figure 1. The curvilinear 
distance along. the intersection of the mid-line of the annulus with 
a meridional plane is denoted by m. The distance normal to the 
mid- line in a meridional plane is labeled by n. The circiim- 
ferential coordinate (j> is considered positive in the counter- 
clockwise direction when viewed down the positive z axis. The 
thickness of the annulus, b, is assumed to be small compared 
to the radius r. Hence, the n component of the velocity vector and 
all variations in the n direction are neglected Cexcept for the 
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regions near the hub and casing). . Accordingly/ the annulus is 
termed a stream surface. The transformation of the governing 
equations to the two-dimensional stream surface coordinate 
system gives the follov/ing results : 

Continuity ; 

* k 

Meridional Momentum; 


9W 3W W. 9W - 

in , /t.t m . (? m W d) . ^2 . ^ 

+ P (W -r— — T ■— T-T -T - - sina - £? r sines - W.sina) 

d t m dm i: d 9 r ^ 


ar. 9 a 3^™ 

= - l£ - 6 i- r„ ( 22. + -M r.A — j. 1 9 


9m 3 9m ^'^e^^m 


r r 94 


■) ] 


* ? ‘35 I T*“'” 


M T 9W, W 
„ "^e .1 4 - m . T 


(7) 


Tangential Momentum : 


aw.' aw, aw. 

pr + p (rW 7 — ^ + W -T-^ 

^ 9t m am 494 


t '^x sina + 2fir W..., sina) 
m 4 m 


.94 3 94 ^9m ~ r 9m“ 


_ _ 


a T W^ aw, 

^4“ ~ sinan+ PgSincs 


W 


(f) * _L 1 3 TT 1 

^ Sina + - -^ W 1 
r 3 9 m 


(8) 
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where 


= relative meridional mean velocity component/ 

=' relative tangential mean velocity component/ 

— 1 <32T 

a = sin ^ = angle between meridional streamline and 
axis of rotation/ fig. 1. 


The system of the differential equations i6) , (7) and (8) , can 
be written alternatively in terms of a stream function i|) and a 
vorticity u. In the present study/ the stream function- 
vorticity formulation is considered for the purpose of simpli- 
fying the handling of the boundary conditions. Among other 
advantages, this formulation offers the ability to express the 
governing equations in the conservation law form. By so doing/ 
the finite difference approximation of the partial differential 
equations will not allow excessive accumulation of the errors 
in the fluxes of the conserved quantities. 

For steady state flow situation/ and for the stream channel 
configuration shown in Figure 1, the continuity equation is 
satisfied by specifying a stream function Ci|)) / with the 
following relations: 


and 


m b pr 3^ 



M 1 3tj) 
b p 3m 


(9) 


m 

where M is the rate of mass flow passing through the stream 
channel. If the pressure terms from equations (7) and (8) are 
eliminated by cross-differentiation and the mean vorticity 
variable w is introduced/ one obtains the vorticity transport 
equarion: 

\ 4-l2. 3p 3W^/2 

3m 3<{) “ 36 ^b 3m 3m 3(J) ” 3^ 3m 

- Ih I s ^ It J " '=1 = » ‘ 
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where 


( 11 ) 


2 2 2 
Vr = VI + KZ 
m <p 

G, = 2J^r (sina) ~ + 2fi(sxna) W. -|x + 2J5r(sina)p 
1 rti 9in © 9 © 

9W 

+ 2J2p (sina) + 25^{sin^a) pW^ + 2firpWj^ 


+ n r sina 


3p 

9<Ji 


( 12 ) 


and <0 is defined by 


to = 


1 

r 



(rW.) 


9 

3 9 


(W )] 

^ m 


(13) 


When the velocities and W^ in equation (13) are expressed xn 
terms of the stream function variable , as defined in equation 
(9), equation (13) would reduce to the stream function equatxon 


to 


i + 1— ( ^ 1 

r ^9m 'bp 9m 9^ 'bpr 9^^ 


(14) 


For a turbomachine rotor it is convenient to express the energy 
equation in terms of the total enthalpy (H) of the gas/ besides 
it velocity components. The total enthalpy for turbulent flow 
is expressed as follows : 

H = h + ^ — + E 

2 ^2 2 

or H = h + |- + + E (15) 


where V = absolute velocxty vector 

E = kinetic energy of turbulence. 

Therefore, the energy equation C3) , when transformed to the two- 
dimensional stream surface coordxnate system, results in 
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where P_ = turbulent Prandtl n\mber = (C u /K) 
r p e 

S^, = turbulent Schmxdt number = (u /r„) 
cn e E 

Tg = exchange coefficient of the kinetic energy of 

turbulence 


(16) 
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and the dissipation function, D, is given by 
■ aw^ 2 , 3W. » 


am" 


B. i — (W r + ) ] } ' 
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The properties of flow passing through the stream surface 
are completely defined by equations (4)‘, (9), (10,), (14) and (16) 

together with the known variations of u^/ K and the given 
boundary conditions. The kinetic viscosity, y , is calculated 
from a two equation model, one expressing the development of 
the turbulent kinetic energy, E, an d.^ the other its dissipation , 
rate, e. These equations are expressed as follows. 
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Turbulent Kinetxc Energy Equation : 

« • 
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ch 
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Dissipation Rate Equatxon : 
3 ,M 3U;. 3 


(M £ 11) 

3m '‘b 3<j)^ 


fM c 11) - 1_ f Q 3- If.) 

3(|> ^b 3m^ 3m 3m' 
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ch 


1 l£) _ r- i- ta 
r 36 ^1 E ^ 




= 0 


(19) 


where and C 2 are constants to be determined from available 
experxmental data. The relation amongst E, s and is given by: 


y 


e 



( 20 ) 


In general, equations (4) through (19) are valid for any 
t\arbomachine geometry or any number of stream surfaces except for the 
two stream surfaces and as shovm in Figure 2, which 

contain the hub and shroud profiles. This exception may be 
attributed to the existence of a large variation in flov? pro- 
perties along the normal, n, to these tv^o surface resulting from 
the presence of the solid boundaries. 

The solution of the above system of equations v^ithin the 
turbomachine passages is carried out nximerically . One can observe 
from equations (10), (14)., (16), (18) and (19), that they are 
coupled elliptic partial differential equations. The highest 
order derivatives in those equations are the second derivatives of 
6/ w, H, E, and e and as such a, H, E and c are the dependent 
variables. From the nature of the problem, none of the terms are 
negligible in the governing equations . The convective terms 
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introduce nonlinearity and also instability if the proper dif- 
ferences are not taken into accomt. Once a solution for these 
variables has been obtained, the velocity distribution can be 
determined from equation C9) . The pressure distribution can then 
be evaluated from either equations CV) or (8). In order to 
solve the elliptic equations by the usual numerical methods, 
it is necessary to define a selected region in the physical domain 
with boundary conditions specified for all the dependent 
variables. 

Description of the Computational Domain in the Physical Space 

The flow region of interest, as shown in Pig. 3, contains the 
blade row and segments of the stream surface, S-, , extending 
upstream and dov/nstream of the row. Due to the circumferential 
periodicity in turbomachine passages, the selected domain need 
to encompass only that fraction of the flow annulus containing 
a single blade to blade passage. The shape and location of the 
periodic boundaries CAB, NH, IH and FG) may be defined 
arbitrarily so long as rheir spacing corresponds to the blade 
pitch 2 tt/Z. The upstream and downstream boundaries (AN, GH) 
are located sufficiently far from the blade so that tangential 
variation along them are ignored- The flov; properties are 
consequently considered to be uniform along the boundary AN and GH. 

Bovmdary Conditions 

In specifying the boundary condirions, two flow situations 
are investigated. Preliminarily, only the case of the laminar 
flow is considered. The turbulent flow case is now in progress. 
Accordingly, in the following specification of boundary conditions, 
no assignment for the boundary values of E and e in equations 
(18) and (19) is needed. Moreover, the flow properties v/ithin the 
turboraachine channels is completely defined through the simul- 
taneous solution of only equations (10) , (14) S'Uid (J.6).. 

Upstream flow boundary AN ; 

It is a common practice in turbomachine flow calculations 
that the magnitude and the direction of the flow velocities , the 
total temperature and the total pressure or density are defined 
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ar the engine inlet. Therefore, along the boundary AN the values 
of ip, 01 , H or their derivatives can be evaluated using the defined 
flow properties. The known magnitude of the inlet relative 
velocity and its direction, as shown in Figxire 3, specify the 
values of 9tji/3m and according to the relations 


3iJ) _ _ bp ^ 34> _ bpr „ 

^ ■ 3 ? - — 


( 21 ) 


Since the inlet stream of the gas is considered to be uniform, 
the absolute value of vorticity oi has to be zero along the 
boundary AN. In a rotating frame of reference, as it is in the 
present case, the relative value of oi is given by the follovring 
expression 

iti = ~ [2Q sinct] 

inler 

The value for the total enghalpy, H, can be defined using the 
specified flow properties. 

The periodic flow boundaries AB, NM and FG, IH ; 

The periodicity condition requires that the magnitude and 
direction of the flow velocity as well as other fluid properties 
be equal at every two corresponding points along AB and MN. 
Similarly, the same condition should apply at every two 
corresponding points along FG and IH. The periodicity condition 
in terms of the dependent variables is satisfied through the 
following conditions. First, equating u, 3a)/3cj), 3ri)/3<!), H, and 
3H/8(|) values on each corresponding points. Second, ensuring 
that the ij) values are differed by unity between the corresponding 
points . 

The blades surfaces boundary (MI and BF) ; 

For the laminar flow case, tvzo boundary conditions over the 
blades surfaces are usually specified. These are rhe non-slip 
condition and the impermeability of the surface in the case of 
blades V7ith no injection. The non-slip condition requires that 
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0 


(23) 


8i|) 

3n 

where n is the normal to the blade surface. On the other hand the 
impermeability condrtxon requires that no component of velocity 
exist in the direction normal to the blade surface. Therefore, 
the blade surfaces are treated as streamlines v/ith the il) values 
specified as zero on the (Ml) surface and vinity on the CBP) 
surface. On either MI or BP surfaces, one has therefore, two 
boundary conditions for * but not boundary condition for w. 

It is a well accepted fact in computational fluid mechanxcs 
to rely on a modified evaluation of equauion (14) to determine 
the boTondary condition for oj. The modxficatxon is introduced 
in an attempt to insure that equation C23) holds, that xs, to 
satisfy the no slip condition. Thxs approach is utilxzed in the 
current study to determine the value of the vorticity, w, over 
the blade surfaces. 

In regard to the thermal boundary conditxons either the 
blade surface temperature is known or the normal derivative 9T/9n 
is specified as zero for the adiabatic wall conditions. In 
either case, equation (15) is used to determine the value of H 
or its derivative that is used as bo^mdary condition for 
equation (16) . 

The downstream boundary GH ; 

Two basic problems arise in the specification of the boundary 
functions for the dependent variables along GH. The first con- 
cerns the behavior of the two dependent variables oi and H. The 
nature of the problem is that, physically, u and H are known only 
if the bounda3:y is located at an arbitrarily large distance from 
the blade surface. In this case the absolute value of vorticity, (o 
is zero, while the value of H is that of the corresponding 
surroundings. The placement of GH ar exceedingly large distances 
from the blade boundary is quite obviously not possible for 
numerical considerations. Therefore, one has to employ some 
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auxiliary condition, usually obtained by experience, to define 
oj and H implicitly. The conditions of zero gradient in the 
meridional direction, i.e. 


3 CO 3H 


(24) 


seems to be appropriate and therefore is employed in the current 
work. 

More important than the specification of the remote boundary 
functions of u and H along GH is the determination of the ij* values 
at the same boundary. The dovmstreara flow velocities, which may 
be used to determine i|» derivatives along GH, and that guarantee 
a unique solution to the problem are not knovm in general apriori. 
Therefore, one has to introduce a supplementary condition, generally 
resulting fcom physical intuition, to define the ^ derivatives. 
Generally, two approaches have been utilized by investigators 
working in the inviscid flow area to deal with this problem. Both 
of them are based on an iteration procedure, through which the Kutta 
condition for tangency of the flow at the blade trailing edge is 
satisfied. This is equivalent to specify the unique solution 
to the problem. Unfortunately, the Kutta condition cannot be 
applied in a realistic manner to solve the present flow 
problem due to its viscous nature. The conservation of angular 
momentum principle is employed as an alternative supplementary 
condition that results in the required unique solution. The 


procedure is as follows: 

Estimated exit flow angles, along GH are used to 

specify the values of derivatives in m direction through the 


following relation 


il. 

9m 


_1 

2tt/z 


tang 


exit 


'exit 


C25) 


The flow field equations are then solved for the boundary function 
of lb given by equation (25). to obtain the velocity and the pressure 
distribution throughout the turbomachine channel. An evaluation 
of the torque developed by the channel is obtained through the 
integration of the difference in pressure and shear forces 
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acting on the blade surfaces. The change in the angular momentuiti 
between the known inlet and the estimated exxt flow conditions 
is calculated. If the value of the predicted torque was not 
equal to the rate of change of the angular momentum then the 
direction of the exit flow velocity is altered. The whole 
procedure is repeated until a satisfactory result is obtained. 


PRELIMINARY RESULTS AND DISCUSSION 

In the present phase of the study, where the case of laminar 
flow is considered, the governing differentxal equations (10) , 

(14) and (16) are reduced to a set of simultaneous algebraic 
finite dxfference equations. These algebraic expressions are 
then recasted into a suitable form for solution by an iterative 
successive substitution technique. A computer program is 
developed for the numerical solution of the resulting algebraxc 
equatxons subjected to the prescrxbed boundary conditions. 

The program is arranged to handle general flow sxtuations 
within turbomachxnery , These may be axial, radxal or mxxed flow 
types. In general, the program requires as an input the shape 
and the confxguration of the stream surface the blades 

geometry, the inlet flow conditions, and the rotational speed 
of the machine. The configuration of the stream surface xs 

usually generated from a merxdxonal flow analysis. The output 
consists of the distribution of the stream fmction, the vorticity 
and the static pressures wxthin the blade passages. The 
variation of merxdional and tangential velocity components from 
blade-to-blade and from the xnlet of the machine to its exit is 
also generated. In cases where blade cooling is consxdered, 
the program has the capability to generate the temperature 
distribution and the local heat transfer coefficients along the 
blade surfaces. 

Six flow cases are investigated usxng the developed program. 
The main purpose was to check the accuracy of the present method 
of analysis in predicring the actual flow behavxor wxthin 
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turbomacbine channels. The accuracy of the method was confirmed 
by a comparison with available experimental data. Five of those 
investigated cases were concerned with inward flow situation 
while the sixth one dealt wiuh an outv/ard flow case. In all 
cases investigated/ the flow was considered to be incompressible 

-5 

and having a constant kinematic viscosity y of 1.789 x 10 

2 ® 

N - Sec/mt . The blade surfaces are assumed to be adiabatic 

with no heat sources or sinks. 

Inward Flow Cases 

These flow cases are those of a radial inflow turbine whose 
rotor consists of 8 radial straight blades. A full description 
of the rotor geometry is given in Fig. 4. The flow patterns 
on a blade-to-blade stream surface are presented over a wide 
range of operating conditions. The blade-to-blade stream sur- 
face is located midway the passage depth of the rotor as 
shown in Pig. 4. Table 1 summarizes the range over which the 
parameters for the numerical solution are varied. The resulrs 
are presented as stream function contour plots/ velocity profiles 
across the rotor passages together with some information con- 
cerning the pressure distribution within these passages. 

The streamline contours for the five inward solution cases 
are shown in Fig. 5, The streamlines are plotted for the region 
between a pair of blades, represented by the heavy thick lines. 

The streamlines are designated by a stream function ratio ’I'/'l' total 
such that the value of a streamline indicates the ratio of flow 
through the passage between the streamline and the pressure 
surface of the blade. Thus, for the given channel configiiration, 
the streamline spacing is indicative of the velocity relative 
to the rotor, with close spacing indicating high velocities and 
wide spacing indicating low velocities . For the operating 
conditions corresponding to case 2, as shown in Pig, 5a, it is ^ 
observed that a recirculating eddy began to form near the pressure 
surface of the blades. As the rotating speed increases, the 
recirculating zone grows much larger, as shown for case 4 in 
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Fig. 5b. The relative velocity near the suction side of the 
blades increases in the later case. This may be attributed to 
the fact that the effective sectional area of the rotor decreases 
with the growth of che recirculating zone- In actual practice/ 
large recirculatory zones are expected to cause higher losses 
in total pressure. Hence/ it is desirable to avoid them through 
an efficient design of the rotor and by selecting a proper 
operating conditions. Prom inspection of Pig, 5 it is clear 
that the size of the recirculating eddy depends upon the relative 
magnitude of the flow rate (M/^p) and the rotor tip speed 
In general/ these zones can be reduced by increasing the mass 
flow rate or decreasing the rotor tip speed. 

The most remarkable feature of the present results is the 
good agreement obtained between the predicted flow behavior and 
the experimental evidence given in reference 2, In all cases 
studied/ the size and the extent to which the recirculating zone 
grows compares favorably well with the experimental data/ as 
indicated in Pigs. 5a/ 5b and 5c. 

Figures 6a/ 6b and 6c show the predicted meridional velocity 
distribution across the blade passages at three different radial 
locations. These locations are selected to correspond to a 
radius of 23.2, 15 and 6.8 cms, respectively. The dotted line 
shown on each figure represents the velocity distribution for 
the inviscid case. From the profile distributions in the 
previous figures, it is evident that near the turbine inlet the 
boundary layers formed over the blade pressure surfaces are 
thicker compared to those formed on the blade suction surfaces. 
One of the unique features of the viscous flow in rotating 
channels is the large variation in the meridional velocity 
profiles as the flow travels downstream in the channel. The 
profile distribution, at stations located away from the turbine 
inlet/ shows that regions of high meridional velocities are 
shifted towards the blades suction surface/ as shown in Figures 
6a, 6b and 6c. ^^hile, regions of high relative meridional 
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velocities are observed to exist near the blades pressure surface 
at subsequent downstream stations. The inviscid flow solutxon 
gives a completely different flow behavior, in this respect, 
when compared with the viscous calculations. Moreover, xn 
some flow situations v/here severe changes take place near the 
rotor tip as xn case 5 in Pig. 6c, the inviscid flow solution 
fails completely even in predicting rhe flow characterxstics , 

All these factors/ in addxtxon to the appearance of the 
reversed flow regions made it clear why the inviscid solutions 
always faxl to produce a realistic prediction of boundary layer 
characteristic parameters in rotating machines when used xn 
conjunction with standard boundary layer interacting techniques. 

Figures 7a, 7b and 7c show the pressure variation between blades 
at different radial locations. These radial locations correspond 
to radius, r = 13.1, 16,3, 19.5 and 22.7 cms. The static 
pressure, p, is plotted in these figures using the nondimens ion al 

P 1 2 

quantity, ^tip' Py' Py mean static 

pressure and density at rotor inlet respectively. As with the 
stream function and velocity results, the present method of 
analysis provides an excellent prediction of the pressure distri- 
bution over a wide range of operating conditions. correlation 
with experimental data of reference 2 as shown in the previous 
figures verify this statement. On the whole, the value of C 
is observed to be larger on the suction side of the blades and 
decreases near the pressure surface. The rate of decrease of 
C values is smaller when the rotational speed is increased. 

Also; in the region lying between the center of the passage and 
the pressure surface, the values become smaller as the 
rorational speed decreases. All these observations imply that 
the static pressure drop from the rotor inleu increases with 
high rotational speed. Near the suction surface, the pressure 
drop increases with the decrease of rotational speed. Such a 
tendency is remarkable, particularly near the rotor inlet. 
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Outward Flow Case 


The capability of the method of analysis to deal with flow 
in diffusion cascades is demonstrated by studying rhe flow 
behavior within the rotor of a radial bladed compressor . The 
compressor rotor under consideration has 12 straight .radial bladesL. 
Other data pertaining to the rotor geometry is shown in Pig, 8. 

The input data for this case is presented also in Table 1. 

A typical distribution of the flow properties on the stream 
siurface located midway (between the shroud and the hub) as 
shown in Pig. 8 are calculated and the corresponding results are 
presented in Figures 9 and 10. As with the inv/ard flow cases, 
the stream function profile of Pig, 9 reflects the appearance 
of the recirculatory eddy. The inviscid flow solution, although 
shows the existence of the negative flow regions as illustrated 
in Pig. 10b, it overestimates the size of the recirculating eddy. 
This overestimation is supported by the existance of large 
negative meridional velocities near the blades suction surface. 

In an actual case, boundary layer phenomena are expected to 
reduce the effective flow area of the passage, thus increasing 
the volume flow rate per unit area through the effective area 
and thereby reducing the size of the eddy. This is exactly 
the same result obtained using the present viscous flow solution 
In conclusion, it may be stated that the present method of 
analysis gives a better representation of the actual flow 
situation within the passage of turbomachine rotors- The 
preservation of the ellipticity of the flow problem is believed 
to be the major reason that results in this better representation. 
The ellipticity is preserved through the consideration of all 
the diffusion terms of the governing equations during the solution 
procedures . The method proves also to be accurate and provides 
invaluable information on the rotor flow characteristics. This 
IS evidenced by the good agreement obtained between the predicted 
results and the available experimental data. The preliminary 
results obtained for the radial bladed passages of turbine and 
compressor had provided a great insight into the flow phenomena 
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and this should encourage the efforts undertaken currently 
to deal with more complicated blade shapes. Such development 
in addition to the inclusion of the turbulence model in the 
method of solution are undoubtably essential for further 
aerodynamic improvement and performance prediction of turbo- 
machines . 
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TABLE 1. 

PARA24ETER 

FOR THE NUMERICAL SOLUTIONS 





Case 

No. 

Designation 

m/t“ /P 
o' o 

n//t“ 

o 

$ at 

Ptin upstream 

(mt/sL) 

AW/ degrees 

3 at 

downstream 
boxmdary 
HG, degrees 

Reynolds No. 

=W ^ . * r . . /v 

min tip 

T 

O 

°K 

P 

° 2 
kg/m"^ 

1 

Inward 

0.0011 

58.5 

26.2 

66.67° 

52.3° 

1.113 

X 

10^ 

288 

10500 

2 

Inward 

0.0011 

87.9 

39.3 

72.27° 

33.7° 

1.175 

X 

10 s 

288 

10540 

3 

Inward 

0.0012 

87.3 

39.3 

74.25° 

37.7° 

1.208 

X 

10^ 

288 

10600 

4 

_ Inward 

0.0012 

116.6 

52.0 

79.77° 

10.5° 

1.219 

X 

10^ 

288 

10670 

5 

Inward 

0.0012 

145.6 

65.4 

81.28° 

6.3° 

1.275 

X 

10^ 

288 

10750 

6 

Oubward 

0.00343 

589.25 

261. 8 

44.5° 

81.2° 

1.370 

X 

10^ 

288 

10330 


Q9 * In all cases studied, the upstream boundary AN of Pig. 3 is located at radius ratio 
r/r^^p of 1.35, while the downstream boundary HG is located at radius ratio 
gg; of 0.25 4. 

I - Meridional Velocity at Inlet. 
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